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The discrete vortex lattices in a ferromagnet/superconductor bilayer are studied when the ferro- 
magnet has periodic stripe domains with an out-of-plane magnetization. The vortices are assumed 
to be situated periodically on chains in the stripe domains. Only up to two vortex chains per do- 
main configurations are considered. When the domain period is fixed, the threshold magnetization 
is calculated at which the transition from the Meissner state to the mixed state occurs. When 
the domain period is not fixed, the equilibrium domain size and vortex positions are calculated, 
depending on the domain's magnetization and the domain wall energy. In equilibrium, the vortices 
in the neighbor domains are half-way shifted, while they are next to each other in the same domain. 
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INTRODUCTION 

In recent years, heterostructures made of type II su- 
perconductors and ferromagnetic pieces have been the 
focus of studies both experimental and theoretical P, 01 ■ 
In these structures, the magnet and the superconductor 
are separated by an oxide layer to avoid the proximity 
effect and spin difffusion that might lead to the supres- 
sion of ferromagnetic (FM) and superconducting(SC) or- 
der parameters. Therefore, the interaction between the 
FM textures and the vortex matter in SC pieces is main- 
tained only by the magnetic fields generated by the SC 
vortices and the magnet. This strong interaction not only 
gives rich physical effects that are not observed in indi- 
vidual subsystems, but also offers new devices that can 
be tuned by weak magnetic fields. 

One of the realizations of such heterostructures is a 
ferromagnetic/superconducting bilayer (FSB). In recent 
years, FSBs have drawn a great deal of attention both 
experimentally and theoretically. On the experimen- 
tal side, nonsymmetric current-voltage characteristics on 
Lao.66Sro.33Mn03/YBa2Cu307_5 were studied It has 
also been reported that stray fields from domain struc- 
tures of an FM film lead to a significant decrease of 
the SC critical temperature Tc in a zero external field, 
whereas is enhanced under an applied field The 
influence of magnetic domain walls on SC critical tem- 
perature was theoretically studied on the basis of the 
Ginzburg-Landau approach [^. On the theoretical side, 
much attention has been paid to the interplay between 
the vortex matter and the FM layer. Sonin studied the 
conditions for the penetration of a vortex near a mag- 
netic domain wall j6|. Helseth et al. investigated the 
pinning of vortices near the magnetic domain walls 
Recently, Laiho et al. investigated the vortex structures 



in the FSB, when the FM layer has domain structure 
with out-of-plane magnetization . Two possible vortex 
structures were shown to occur: vortices with alternating 
directions corresponding to the direction of the magne- 
tization in FM domain and vortex semiloops connecting 
the adjacent FM domains. 

Earlier Lyuksyutov and Pokrovsky 0,0] noticed that, 
in a bilayer consisting of homogeneous SC and FM films 
with the magnetization normal to the plane, SC vortices 
occur spontaneously in the ground state, even though the 
magnet does not generate a magnetic field in the SC film. 

In previous work, we presented a theory of such vortex- 
generation instability and the resulting vortex structures 
[ll|. We showed that, due to this instability, domains 
with alternating magnetization and vortex directions oc- 
cur in a FSB. In that study, the domain structures were 
treated in the continuum regime in which the domain size 
was much larger than the effective penetration depth, 
A — X^/dsc, where the London penetration depth A is 
much larger than the thickness of a SC film dsc [l^ . 
Under the continuum aproximation, the energy of stripe 
phase was found to be minimal. The equilibrium domain 
size and the e quil ibrium energy for the stripe structure 
were found as 
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where m — m — e^/^o, £v — {(f^o / ^(^'^'^ ^) ^"^{^ / is the 
self-energy of a vortex, Sdw is the domain wall energy per 
domain wall length, A is the domain's area and C ^ 0.577 
is the Euler-Mascheroni constant. If Sdw < 4to^, the 
continuum approximation becomes invalid, since L^q be- 
comes on the order of or less than A (see Eq.Q). How- 
ever, it can be recovered by considering the discrete lat- 
tice of vortices instead. 
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In this paper, the discrete lattice of vortices in the 
stripe domains is studied through a method that works 
in both continuum and the discrete regimes. The method 
we use here is based on London-Maxwell equations and 
was developed elsewhere 0. Earlier, its extention to 
periodic systems for the case of square magnetic dot ar- 
rays on a SC film was used In this work, it is 
adapted to the discrete vortex lattices in SC/FM bilay- 
ers. In doing so, we assumed that vortices and antivor- 
tices sit periodically on chains in the alternating domains 
of magnetization and vorticity. Recently, Karapetrov et 
al. observed vortex chains in mesoscopic superconductor- 
normal metal hetrostructures by means of scanning tun- 
neling spectroscopy techniques Vortex chains are 
also observed in anisotropic high Tc superconductors (see 
review We first studied the vortex lattices in the 
FSB when the FM domain size is fixed. In this case, we 
investigated when the vortices appear spontaneously in 
the FSB, depending on the domain size and the magneti- 
zation strength of the domains. We found that one chain 
per domain configuration appears first, followed by two 
chains per domain state upon further increase of magne- 
tization at the fixed domain period. Next, we considered 
the case in which the domain size is not fixed. For this 
case, the following problems were studied; i) when the 
domain structures appear spontaneously; ii) how the vor- 
tices and the antivortices are positioned on the chains; iii) 
how the equilibrium domain size changes, depending on 
the magnetization and the magnetic domain wall energy 
in the presence of the vortices. In order to solve these 
problems, we first proposed five different configurations 
of the vortex and the antivortex chains, in which at most 
two chains per stripe are considered. Next, we calcu- 
lated their equilibrium energies by means of numerical 
methods and found the most favorable case among them 
by minimizing their energies with respect to domain size 
and vortex positions. Our calculations showed that, in 
equilibrium structure, vortex chains are half shifted in the 
adjacent domain, while they are next to each other in the 
same domains. The comparison of equilibrium energies 
of cases with one and two chains per domains shows that, 
at lower values of magnetization and domain wall energy, 
the case with two chains is energetically favorable. Addi- 
tionally, the single-chain case does not win over the ones 
with two chains per domain in the equilibrium domain 
structures. 

The outline of this paper is as follows: in the follow- 
ing section, the method for the discrete case and its ap- 
plication to configurations with single-vortex chain and 
double- vortex chains per domain is introduced. The third 
section is devoted to our results for vortex chain states in 
the domains with a fixed period. In the fourth section, 
we present our results on the proposed five configurations 
of single-chain and double-chain cases when they are in 
the equilbrium domain structures. The last section con- 
sists of the conclusions and discussion. In the appendix. 



the details of the methods and mathematical tricks in the 
series calculations are given. 

METHOD 

In the continuum approximation, we found that the 
vortex density increases at closer distances to the mag- 
netic domain walls. Based on this fact and the symme- 
try of the stripe domain structure, it is reasonable to 
consider that the vortices and antivortices form periodic 
structures on straight chains along the y direction. Even 
though it is not clear how many chains are associated 
with each domain, we can still make progress toward un- 
derstanding discrete vortex lattices. To this end, five 
stripe domain configurations in which vortices are situ- 
ated periodically on chains are proposed. From this point 
on, the configurations with N vortex chains per stripe 
domain are labeled as N state. 

In two of the proposed cases, there is one chain per 
stripe [N ~ 1 states), located in the middle of the do- 
main. In this case, two configurations of vortex lattice 
are possible. First, the vortices and the antivortices in a 
neighboring domains are alongside one another (see Fig. 
1(a) I; second, they are shifted by half period b/2 along 
the y direction, where h is the distance between two near- 
est vortices on the chain (see Fig. 1(b) I. 

In the other three candidate lattice structures, there 
are two chains per stripe domain N = 2 states, at a 
distance a from the magnetic domain walls. The possible 
situations in the two-chain cases are as follows. First, 
chains in the same domains are shifted by a half period 
6/2 along the y direction (see Fig. 



1(c) I. However, the 



neighbor vortices and antivortices are next to each other. 
Second, they are just shifted by a half period b/2 (see Fig. 
1(e) I. Third, the vortices and the antivortices are simply 
side by side on the chains (see Fig. 1(f) I. 

Our next step was to write the energies of these five 
candidates. To this end, we used the en ergy equations 
for periodic systems developed elsewhere |l4j . 
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In these equations, the vortex configurations differ by 
their form factors. We can obtain them from Fq = 
J2r nie^^'^\ where the G's are the reciprocal vectors 
of the periodic structures, the are the positions of the 
vortex centers, and rii are the charge of the vortex. In 
our proposed models, G = ((2r-f 1)?-, 2s?) and rii = ±1. 
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TABLE I: The form factors of vortices in the proposed con- 
figurations. 
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(a) 1st Case (b) 2nd Case 



Fig. 1(d) I. Since information about the vortex lattice 
is carried only by the form factors, there is no need to 
consider the above-mentioned case separately. 

In our calculations, the divergent part of the series 
must be extracted carefully. We show below a detailed 
analysis of the series equations for each candidate. We 
start with the self-interaction energy of the magnetic 
layer Umrm since it is the same for each configuration. 
For the periodic structures, it is given by Eq. Q. Direct 



substitution of the Fourier coefficient of the stripe phase 
nT-zG = ^(2r+i) ^"^^^ ^'i- © gives the self-interaction of 
the magnetic layer per unit cell as 
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where ip^^\x) is the polygamma function of zeroth order 
|l7|. In our numerical calculations, we write the loga- 
rithmic term in Eq. © as ln(A//) -I- ln(L/A) and then 
incorporate the — 4m^ln(A/^) term in the renormalized 
Another energy term with a divergent series is the 
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vortex energy, in general given by Eq. J^Jl. The loga- 
rithmic divergence in this term stems from the vortex 
self-energies. We first split Eq. © into two parts as 
follows: 
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(e) 4th Case 
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FIG. 1: Proposed configurations for A'^ = 1 and N — 2 states. 



TablelH gives the form factors of each configuration in the 
order they are described above. 

Note that, in Table HJ the form factor for the third con- 
figuration also belongs to the case in which the vortex 
and the antivortex chains are shifted by half period only 
in the neighbour domains, not in the same domain (see 
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Note that the area of the unit cell is 2Lb. The first term 
of the series above contributes to the self-energies of the 
vortices; whereas, the second term is the vortex-vortex 
energy and will be left in the series form. For each form 
factor in Tabled the series in the first term can be trans- 
formed to the form of 

J2T=-oo l/((2r+l)^a;^ + s^), where x is constant, 

and depends on the form factor. A detailed analysis of 
such a series is given in the appendix. 

The next step is to find the vortex energy and the 
interaction energy of the magnetization and vortices 
for each configuration. In the calculation of Umvi we 
take the Fourier coefficient of the magnetization to be 
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(2r'+i) ^i^v)- ^^'^ i^^&i the stripe is infinite along 
the y direction results in the additional term 27r(5(Gy). 
However, it does not play any role in the calculation of 
Umm- For numerical analysis, these energies must be ex- 
pressed in terms of dimensionless parameters. To this 
end, we define dimensionless variables A = A/L ,b = b/L 



In the third configuration, as in the second configura- 
tion, only even values of s contribute to the energy. In 
the first two configurations, the square of their form fac- 
tors enter the vortex energy as a constant. However, in 
this case, the square of the sine function appears. In the 



and Sdw = ^dw-^/^o- The total energy t/ is measured in 
units of Eq/A^. In addition, we introduce the dimension- 
less magnetic energy as Umm = Umm/i^o/J^)- In terms 
of these parameters, the energy of the first configuration 
reads 
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appendix, the calculation of the series in the presence of 
such functions is shown. Introducing the dimensionless 
parameter a = a/L^ the energy functional of the third 
configuration becomes 



where. 



^ ^ coth((2r + l)^)-l 
Jv 2^ 2r + l 

oo ^ 
fvv ^ ^ 



r,s= — oo 

oo 



((2r + ly + + 27rA^(2r + l)2 + ^)) 



" ^(2r- + l)(l + 27rA(2r + l))' 



I 

The form factor for the second configuration survives of the second configuration is found to be 
only for even values of s. Then, the dimensionless energy 
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where, 



+ Umm + £dwK (12) 
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In the fourth configuration, the square of the form fac- of s give different contributions. Then, we can calcu- 
tor is: late the vortex energy for even s and odd s separately. 

|Fg|^ = 2 — 2(— l)*cos((2r + l)7ra). Even and odd values Employing similar techniques, we find 



26 V Ve'-^AC/ ^i" ^0 / 



where f^l = f^l and. 
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The form factor for the fifth case resembles that of 
the third case with an exception. That is, in the third 
case, only even values of s are taken into account, while 



all integers contribute to the sum over s in the fifth case. 
Keeping this in mind, obtaining the dimensionless energy 
for the last case is straightforward: 



U''' = |(ln(-ii^)-ln(cotM)+4/(^)(A,a)-i-/(^)(A,a^ 

+ Umm + Sdw^, (16) 



where fml = fml and. 
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COMPUTATIONAL METHODOLOGY 

The series in /i*' converges very fast when r^ax > 200, 
while the series in fiv converges rather slowly. When 
Tmax > 4000 and s„iax > 4000, the results do not change 
up to the 6th decimal point in the energy, where i labels 
the particular domain configuration. To make sure of this 
accuracy in the calculations, we take r^^ax = 600 for 0\ 
Tmax = 5000 for fly}v, and rmax = 5000 and Smax = 5000 
for f'ir} in Eqs.ljHmniinilllllSl)- In addition, the respec- 
tive error deviations for A/L,6/i and a/L in numerical 
calculations arc ±0.005, ±0.0005 and ±0.0125. In the nu- 
merical computations of Eqs.||H| El El El I16|l . we take 
ln(4A/(e^0) = 5.57. 



THE VORTEX CHAIN STATES IN DOMAINS 
WITH FIXED PERIOD 

Here, the conditions for the transition from the Meiss- 
ner state to the mixed state when the vortices first 
spontaneously appear, are determined. As described in 
the second section, the total energy of the bilayer con- 
sists of the self- vortex energies , vortex- vortex interac- 
tion Uy^j , vortex-magnetization interaction Umv , the self- 
interaction of magnetic domains Umm and the domain 
wall energy Udw The last two terms are necessary only to 
determine the equilibrium domain size. Therefore, they 
are irrelevant when the domain period is fixed, and they 
will be ignored in the further calculations to determine 
when the Meissner state-mixed state transition occurs. 
Now, we call the total energy of interest, the effective en- 
ergy to avoid any confusion in the rest of the paper. The 
effective energy has the form U = Uy + Uyy + The 
detailed versions of these terms for each proposed config- 
uration were described in the second section. Next, we 
calculate the equilibrium effective energies of each con- 
figuration for fixed values of m^o/eo and A/L. In doing 
so, the effective energies are minimized with respect to 
the vortex positions. The necessary condition for the 
vortices to appear spontaneously is C/e// < 0. Our nu- 
merical calculations show that the configurations with a 
single vortex chain per domain (N = 1 state ) appear 
first in the mixed state. Note that there are two configu- 
rations for the = 1 state. In our calculations, it turns 
out that they are indistinguishable when the transition 
occurs. Next, the values of m0o/eo that make C^j/J^ zero 
for various values of the domain period, are calculated. 
Thus, the curve that separates the Meissner state and 
= 1 state is obtained (see the bottom curve in Fig. 
121). The fit of our numerical data yields an equation for 
this curve as: m^o/eo = 2.75(A/L)'^-^^. Earlier, Laiho 
et al. calculated critical magnetizations for two cases: 
penetration of a vortex from a domain center and pen- 
etration of a vortex semiloop from a domain wall in a 
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FIG. 2: Phase diagram showing transition from the Meissner 
state to the mixed state, depending on m(f)o/eo and A/L. The 
region below the TV = 1 curve corresponds to the Meissner 
state, while N = 1 state exists in the area between the A = 1 
and N = 2 curves. The N = 2 state becomes energetically 
favorable above the N — 2 curve. 



FSB, in which both FM and SC layers are thick. For the 
former case, they found that the critical magnetization 
decreases with the increase of domain period. This result 
is in full agreement with ours. When vortices appear first 
in FSB, the inter-vortex distance b grows. With further 
increase of the magnetization, they get close. Note that 
our results here do not give any information regarding 
which configuration of A = 1 becomes favored in the 
mixed state. 

We also studied the transition from the A = 1 state 
to the N = 2 state. The necessary conditions for this 

transition are U^X < ^efT ^^'^ < ^- 

culations show that the third configuration wins over the 
fourth and fifth configurations for the N = 2 state. Note 
that the effective energy of the third configuration corre- 
sponds to two cases that have the same structure factor. 
Which one is likely to win over is discussed in the next 
section. Following steps similar to those described above, 
we obtained the curve m^o/eo — 4.72(A/L)° ''^ (see the 
top curve in Fig. Furthermore, when the N — 2 

state wins over the A = 1 state, b and a are on the or- 
der of a few A. The increasing magnetization leads to a 
significant decrease of these parameters. 
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THE VORTEX CHAIN STATES IN THE 
EQUILIBRIUM DOMAIN STRUCTURES 

Here, we first investigated when the proposed cases 
become energetically favorable. To this end, we checked 
when the equilibrium energies of the cases first become 
negative. To do so, the energies of five proposed cases 
Eas. H8linil2ll4ll(T|l for different values of m(j)o/eo and idw 
were calculated by minimizing their energies with respect 
to the domain size and the vortex positions. From this 
procedure, one can determine which case's equilibrium 
energy becomes negative first and then calculate the cor- 
responding values of mipo/eQ and idw Our calculations 
show that the equilibrium energy of the third case associ- 
ated with the N = 2 state turns to negative first. When 
the N = 2 state first appears, the domain size and in- 
ter vortex distance on the same chain are on the order 
of a few tens of A. This suggests that, in equilibrium, 
vortex chain states appear near the domain walls. It is 
expected that vortex chain states proliferate with further 
increase of the magnetization. This result is consistent 



with our previous work in which the vortex density in 
the continuum aproximation increases substantially near 
the magnetic domain walls. In Fig. 01 the phase dia- 
gram for the N = 2 state is depicted. The dots in the 
figure are obtained from our simualtions. The curve is 
the fit to the simulation data. The fit gives the curve 
m0o/eo = 1.37e°i;^'^. Below the curve, there is not any 
stable configuration, while, above the curve, the third 
configuration corresponding to TV = 2 state exists. 

Nonetheless, this information is not enough for us to 
understand the equilibrium structure, since corre- 
sponds to two different configurations with the same 
structure factor. At this point, further analysis is need- 
edto determine which configuration is more likely. To 
this end, we calculated the minimal energy of each case 
in both continuum and discrete regimes, depending on 
m0o/eo at fixed idw The equilibrium energies for these 
cases in discrete and continuum regimes are given in Ta- 
ble^ Note that small and large values of m(/)o/eo for 
each Edw correspond to continuum and discrete regimes, 
respectively. 



TABLE IL Equilibrum energies for proposed configurations. Two columns on the left are input. 





m(j)o/£v 


f/i 


r/2 


t/3 


U4 


f/5 


0.01 


5 


-2.58454176 


-2.58455668 


-3.36407195 


-3.36404625 


-3.36404615 


0.01 


20 


-65.98296440 


-65.98296500 


-89.20105311 


-89.20030961 


-89.20030943 


0.1 


5 


-2.54211623 


-2.54211637 


-3.33057991 


-3.33054828 


-3.32949063 


0.1 


20 


-65.87136440 


-65.87136500 


-89.10475311 


-89.10310961 


-89.10310943 


1 


5 


-2.16635495 


-2.16637343 


-3.01972350 


-3.01972355 


-3.01972347 


1 


20 


-64.78187766 


-64.781878054 


-88.15826391 


-88.15826394 


-88.15826367 


10 
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-0.35644861 


-0.35667498 


-1.21258278 


-1.21232397 


-1.21232396 


10 


25 


-95.00241100 


-95.00247780 


-134.28090345 


-134.27589780 


-134.27589769 



In our numerical calculations, we find that all proposed 
configurations are stable in both the discrete and contin- 
uum regimes, indicating that our method works well in 
both regimes. As in the previous section, the third case 
associated with the N — 2 state wins over the other two 
chains-per-domain configurations. As said before, the 
energy of the third case corresponds to two configura- 
tions. Which configuration is more likely to appear in 
equilibrium can be figured out from simple physical con- 
siderations. Namely, in FSB, the equilibrium structure 
is determined by the competition between vortex-vortex 
and vortex-magnetization interactions. The former fa- 
vors vortices and antivortices in the neighbor domains 
to line up in a transverse direction (perpendicular to the 
magnetic domain wall), whereas the latter prefers vor- 
tices and antivortices to be shifted so that gain in energy 



is maximized. When vortices are next to each other on 
either side of the magnetic domain wall, the magnetic 
fields they produce cancel out each other. From the nu- 
merical results, it is obvious that vortex-magnctization 
interaction wins the competition and results in half-way 
shifting of vortices, if one compares the energies of the 
1st and 2nd cases. Then, vortex-magnetization interac- 
tion is the dominant factor. By the same token, one 
can understand what is going on in double-vortex chain 
configurations. For instance, in the fifth configuration, 
energy gain due to vortex-magnetization interaction is 
diminished, since all the vortices are side by side. This 
explains why the equilibrium energy of the fifth configu- 
ration is higher than those of the third and fourth cases. 
In the fourth configuration, the vortices and antivortices 
in the neighbor domains are shifted half-way, so that this 
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FIG. 3: The phase diagram for the third case N — 2 state. 
The N = 2 state becomes stable in the region above the curve. 



configuration must be preferred over the one in which 
they sit side by side in the neighbor domains according 
to the above arguments. However, an alternative config- 
uration for the third case has two chains shifted half-way 
in the neighbor domains instead of one chain as in the 
fourth case. Therefore, one might expect the gain to be 
even more than that in the fourth configuration. An- 
other interesting result is that the system does not favor 
the iV = 1 state at all. Actually, this does not surprise us, 
since, in the continuum approximation, we found that the 
vortex density increases near the magnetic domain walls. 
This fact already suggests that the system favors vortex 
chains being near the magnetic domain walls rather than 
a single chain in the middle of the domain. 

In numerical calculations, equilibrium domain size 
L/A, vortex-vortex distances on the same chain 6/A and 
vortex-magnetic domain wall distances a/ A are also cal- 
culated. Results for the third configuration at various 
values of m(j)o/eo and Sdw are depicted in Figs. ^ 
El Results for other configurations are not shown here, 
since they look quite similar. At fixed Sdw, the further 
increase of riKjjo/^o shrinks the domain width, while a 
higher domain wall energy favors a larger domain width 
at fixed mcjjo/so, as in usual ferromagnets. That is to 
say, ferromagnet favors narrower domains to minimize 
the demagnetization energy, whereas domain wall energy 
makes domains wider. The competition between these 
two energies determines the domain size. Here, the pa- 
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FIG. 4: L/A versus m0o/eo for the third configuration. 
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FIG. 5: 6/A versus m^o/eo for the third configuration. 



rameter m^o/^o plays the role of demagnetization en- 
ergy. Domain wall energy does not affect the distance 
between the vortices located on the same chain. However, 
at larger values of m(f)o /sq, the vortices on the same chain 
get closer. This implies that the unit cell area shrinks, 
and consequently vortex density per area increases. 
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FIG. 6; a/A versus m<j)o/eo for the third configuration. 
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CONCLUSIONS 

In this article, we reported our results on the lattices of 
discrete vortices in stripe domains in FSB using a method 
based on Maxwell-London equations. If Sdw < 4m^, the 
continuum approximation becomes invalid. Instead, we 
considered the discrete lattice of vortices in which the 
vortices were considered to be situated on chains directed 
along the stripes. We analyzed the vortex configurations 
up to two vortex chains in two cases: first, we took the 
domain period to be fixed, second, the domain period was 
not fixed. In the former case, we calculated the threshold 
magnetization at which the transition from the Meissner 
and mixed state occurred. We showed that configura- 
tions corresponding to the = I state appeared first in 
the mixed state. With a further increase of magnetiza- 
tion, the N = 2 state becomes more favorable. The crit- 
ical magnetization for transition from the A'^ = 1 to the 
N = 2 state was also calculated. In the second case, de- 
pending on the magnetization and the magnetic domain 
wall energy, the equilibrium energy, the positions of the 
vortices and the equilbrium domain size were calculated. 
According to our calculations, in equilibrium, vortices on 
the either side of the magnetic domain walls are not side 
by side on the chains; instead, they are shifted by a half 
period along the stripe, while they are side by side in the 
same domain. The threshold magnetization as a function 
of domain wall energy for this case was calculated. 

In numerical calculations, we also found that the vor- 
tex lattice is stable for Sdw > 4m^. At this point, the 
domain size is noticably larger than the effective pene- 
tration depth A, so the continuum approximation is valid. 
Therefore, we expect that the domain nucleation starts 
in the continuum regime. This problem is left for the fu- 
ture research. At constant idw, with increasing rrKpo/sy, 
the equilibrium size of the domain decreases. In addition, 
the vortices on the chain get closer to each other. These 
results agree with those obtained in the continuum ap- 
proximation. As ediu/4m^ increases, we expect that new 
vortex chains develop within the domains. We leave the 
detailed analysis of this problem to another publication. 
In concluding, we believe that confined geometries, such 
as domains in FSB, lead to novel vortex structures, which 
might be experimentally investigated by means of scan- 
ning tunnelling spectroscopy. 
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CALCULATIONS OF SERIES 

In this appendix, the detailed analysis of the series is 
given. First, the series in the energy calculations of the 
periodic systems are analyzed; second, the detailed cal- 
culation of the vortex density is shown. The series we 
encounter in the energy calculations fall into two cate- 
gories. In the first category, we sum over one variable. 
The series in this category are in the form of X)r=i°^ 1/r . 
Employing the Euler-Maclaurin summation formula |l8| , 
the summation is found with logarithmic accuracy as 



In 7 



(A.18) 



where C ^ 0.577 is the Euler-Mascheroni constant. If 
the summation is performed over only odd integers, we 
can still transform our series to Eq. lf08)l . Namely, 



^ 2r + 1 ~ ^ r ^ 2 ^ r' 



ln(2w + l)-fC-ln(^)-|, 
1 



(lnr,„„^ + C + 21n2). 



(A.19) 



The other double series of interest here is in the form 



of 



S — OO 



(A.20) 



where x is an arbitrary constant. Although Eq. ljA.20|) 
is logarithmically divergent, the sum over one of the 
variables can be done easily. To this end, we perform 
the sum over s first. In doing so, Eq. HA.20|I becomes 
(27r/x) coth(7rxr)/r This series is logarithmi- 

cally divergent. In order to get the logarithmic term , we 
add and subtract 1/r. Using the result in Eq. (|A.18|I . 
finally we get 
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2-K 



E°° coth(7rxr) — 1 
+1^'^" 



C 



(A.21) 



Employing the same techniques, we give the results of 
the different versions of Eq. (fOol below: 
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r=oo s— oo 

E E 

r— — oo s— — CX2 



1 



27r 



a;2(2r + 1)2 + 



.r=0 



coth((2r + l)nx) - 1 
2r + 1 



(A.22) 



r— OO s— oo 
r— — OO s— — oo 



1 



x2(2r + l)2 + (2s + l)2 



^ tanh((2r + l)i^)-l Inw C 
^ 2r + l 2 2 

,r=0 



r 



(A.23) 



In Eas. (|I^ and lfOS|) . we use E^ol/(y^ + (^^ + similar way, using sin^((2r+l)2/) = (l-cos(2(2r+l)?/))/2 
1)2) = 7rtanh(7ry/2)/(4?/). In the presence of sin^((2r + or cos2((2r+ l)y) = (l + cos(2(2r+l)y))/2. For example, 
1)?;) or cos2((2r + 1)?/), the series can be calculated in a 



r=co s— oo 

E E 

r— — oo s — — oo 



sin2((2r + l)y) 27r 



E 



sin2((2r + l)y){coth{{2r + l)nx) - 1) 
2r + 1 



(a;2(2r + 1)2 + s^) x 

, Inrmax In|cot(y/2)| C 



(A.24) 



r— OO s— oo 

E E 

r— ~oo s — — oo 



cos2((2r + 1)?;) _ 27r 
(a;2(2r + 1)2 + s2) " ~ 
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We use 



>^ sin2((2r + l)?/)(coth((2r + l)i:x) - 1) 
2r + 1 

r— 

^ ln|cot(y/2)| , C 



E 

r=0 



cos((2r + 1)6*) _ ln|cot(6'/2)| 
2r + 1 ~ 2 ■ 



(A.25) 



(A.26) 



r 
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